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Abstract 

We report on calculations of the translational and rotational short-time self- 
diffusion coefficients D* and for suspensions of charge-stabilized colloidal spheres. 
These diffusion coefficients are affected by electrostatic forces and many-body hydro- 
dynamic interactions (HI). Our computations account for both two-body and three- 
body HI. For strongly charged particles, we predict interesting nonlinear scaling re- 
lations DI (X 1 — at4>^/^ and -D^ oc \ — a^cfp' depending on volume fraction <j), with 
essentially charge- independent parameters at and ar ■ These scaling relations are strik- 
ingly different from the corresponding results for hard spheres. Our numerical results 
can be explained using a model of effective hard spheres. Moreover, we perceptibly 
improve the known result for of hard sphere suspensions. 

Keywords: Self-diffusion, Hydrodynamic Interaction, Charge-stabilized Colloidal 
Suspensions 

PACS: 82.70.Dd, 83.10.Pp 

Self-diffusion of spherical colloidal particles has been studied experimentally over a wide 
range of time scales by means of various scattering techniques, in particular by polarized 
and depolarized dynamic light scattering (DLS). At short times on the scale of DLS, the 
particles have only moved a small fraction of their diameter a, and the particle motion is 
determined by solvent-mediated many-body hydrodynamic interactions (HI) weighted by 
the equilibrium microstructure. The latter is determined by direct potential forces arising, 
e.g., for hard-sphere particles from the steric repulsion between the particles, and, in the case 
of charge-stabilized particles, from the electrostatic repulsion of overlapping double layers 
0. The configuration-averaged effect of HI gives rise to values of the translational and 
rotational diffusion coefficients D\ and D'l that are smaller than their respective Stokesian 
values at infinite dilution, i.e. Dq = ksT/ {Gnrja) and Dq = ksT/ {Snrja'^). Here, a is the 
particle radius, and t] is the shear viscosity of the suspending fluid. 
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The properties of hard spheres are in various respects easier to describe quantitatively 
than those of charge-stabihzed particles. As a consequence, there are many experimental 

1^, Q and theoretical |, |^ results available on the short-time self-diffusion coefficients 
of hard spheres. With regard to the computation of the first and second virial coefficients 
of D\ and in an expansion in terms of the volume fraction 0, the currently established 
results for the normalized diffusion coefficients Hi and HI are given by 0, |^ 

Hi = = l-1.831(/) + 0.88</)2 + O((/)3), (1) 

Dq 

HI = = 1 - 0.6300-0.6702 + 0(0=^). (2) 

Dq 

The possibility to expand Hi and iJJ in powers of arises from the fact that hard-sphere 
suspensions at small are dilute both with respect to HI and to the microstructure. 

While the short-time dynamics of hard spheres is well understood, far less is known thus 
far about charge-stabilized suspensions. The purpose of this letter is to show that there 
are striking differences in the 0-dependence of Hi and H^ between charged and uncharged 
suspensions, and also to provide quantitative predictions. These unexpected differences are 
most pronounced for deionized, i.e. for salt-free suspensions of charged particles. For such 
systems, our numerical results for if* and i7J are well represented by the parametric form 
1 + p0°, where a is an exponent larger than one. Due to the strong direct interparticle 
interactions, deionized suspensions especially exhibit pronounced spatial correlations even 
for very small 0, say < 10~^, so that contrary to hard spheres these systems are dilute 
only with regard to HI. The corresponding radial distribution function (rdf) g{r) has a well 
developed first maximum and it exhibits a so-called correlation hole, i.e. g{r) is essentially 
zero up to a well-defined nearest-neighbor separation larger than cr . In comparison, the 
rdf of hard spheres is nearly equal to a unit step function B(r — a) for < 0.05, and an 
analytical expression for g{r) of hard spheres is known up to first order in [|]. Therefore, 
the calculation of if* and H^ at small is more demanding for charged suspensions, because 
it is necessary to use static distribution functions generated by integral equation methods 
or computer simulations. 

We base our calculations of Hi for charge-stabilized suspensions on the general expres- 
sion Hi = (TrD"^(r^))/(3i5Q) as derived from the generalized Smoluchowski equation 0. 
The corresponding expression for H^ is obtained by replacing the superscipt t by r. The 
hydrodynamic diffusivity tensor D**L(r^) (D'J^(r^)) relates the force (torque) exerted by 
the solvent on an arbitrary particle 1 with its translational (rotational) velocity. Tr D**^ 
denotes the trace of Hfi, and the factor 1/3 accounts for spatial isotropy. Due to the many- 
body character of HI, both tensors depend on the instantaneous A^-particle configuration 
r-^ = (ri, . . . ,rjv), and in principle the full A^-particle distribution function is needed to 
perform the equilibrium ensemble average (...). Thus, it is not possible to perform an exact 
calculation of if* and ifj, that is valid for all particle concentrations. For small 0, however, 
when the mean particle distance gets sufficiently large, a good approximation for these 
quantities is obtained by considering only two-body and, to leading order, three-body con- 
tributions to the HI. For this reason, we use a rooted cluster expansion for the calculation 



of Hi P, |T0[ , leading to the following series expansion of H, 



Hl = l + Hl<P + Hl^<P' + ... , (3) 
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which we truncate after the third term. Here, Hl^ is given in terms of integrals 

Hi, = - dr r'g{r)67rr^a (a\\{r) + 2/3f,(r)) , (4) 
involving g{r) and scalar two-body mobility functions «fi(r) and f3l\{r), whose expansions 



in powers of (a/r) are known, in principle, up to arbitrary order ||TT] , |T2|. In our calculations, 
we include contributions to afj^ and l3l\ up to (9(r~20). The coefficient if*2 is far more 
difficult to calculate since it involves three-body HI. By considering the leading term in the 
far- field expansion of the three-body part of D*\ , i?*2 is approximated by the threefold 
integral 

225 /"i /"i /"i 

HI2 = dti2 dhs d^ g'^'\tu,tn,0ftiti2,ti3,0, (5) 

o4 JO JO J-i 

/t(ti2,ti3,0 = ^^|ii44- 2 (4 + 4) -10^^12^13(4 + 4) 

+ {lotj.ti, + 6 (4 + 4)) - eetutis (t?3 + t'u) + s^'ti.ti,^, 

with /i(ti2, ^13, = ^12 + '^13 — 2^ti2ti3. This integral involves the static triplet correlation 
function g^^^ expressed in terms of ti2 = 2a/ri2, ^13 = 2a/ri3, and ^ = ■ ri3/(?"i2'^i3)5 
where = — is the relative vector between the particles i and j, and is its 
magnitude. 

A similar analysis is used by us for calculating ifj, leading to expressions for the co- 
efficients and HI2, which appear in a series similar to Eq. (^), and which involve 
now rotational two-body and three-body mobility functions. For conciseness, we will not 
quote here the expressions for H^, and HI2, since these are given in Ref. 0. Once again, 
we account for terms up to 0{r~'^'^) in the far-field expansion for the two-body mobility 
functions, and for the leading three-body part of D^^ 

For charge-stabilized suspensions, it is only necessary to account for the first few terms 
in the expansion of the two-body mobility functions, since the integrals in Eq. converge 
rapidly because g{r) is practically zero at small r 0. On the other hand, many terms are 
needed for hard spheres to accurately obtain the first virial coefficients as depicted in Eqs. 
(|l]) and (^. Notice that the second virial coefficients for hard spheres are made up of two 
contributions. The first one is due to Hi, and H^,, with g{r) expanded up to first order in 
0, whereas the second one arises from three-particle HI as embodied in HI2 and ifjg B 0- 
The second virial coefficient of in Eq. was obtained by essentially accounting for 
all two-body contributions in H^,, and also for the leading three-body contribution 0]. On 
the other hand, only two-body terms up to (9(r~^) plus the leading three-body term were 
used so far in calculating the second virial coefficient of Hi, as given by the value 0.88 in 
Eq- (0) 0- By considering terms up to (9(r~20) calculating Hi,, we obtain an improved 
value of —1.096 for the two-body part of the second virial coefficient. Together with the 
three-body contribution HI2 = 1.81, which is obtained by Monte-Carlo integration of Eq. 
(H), we get the improved result 

i/* = 1-1.8310 + 0.710^ (6) 

This result is in better agreement with experimental data ^ for hard-sphere suspensions 
than Eq. (^. The experimental data in Ref. 0] especially agree almost perfectly with Eq. 



However, for charge-stabilized suspensions, it is not possible to use low-order virial 
expressions of the static distribution functions. We use instead results for g{r), obtained 
from the rescaled mean spherical approximation (RMSA), as applied to the one-component 
macroion fluid model (OCM) of charge-stabilized suspensions [0. In the OCM, the effective 
pair potential u{r) acting between two particles consists of a hard-core part with diameter 
a, and of a screened Coulomb potential (3u{r) = Kaexp [— ^(r — cr)] /r for r > a. Here, 
K = (Ls/cr)Z^(H-/«cr/2)^^, Lb = e'^/^eksT), e is the dielectric constant of the solvent, and 
Z is the effective charge of a particle in units of the elementary charge e. k is given by the 
Debye-Hiickel relation = Lb ['24:\Z\(j)/a^ + Simg], where Ug is the number density of added 
1-1-electrolyte, and the counterions are assumed to be monovalent 0. For computing 
and HI2, g^^\r,r') is needed as static input. To this end, we use for simplicity Kirkwood's 
superposition approximation for ^''•^''(r, r'), with the rdf calculated in RMSA. The threefold 
integrals are calculated using a Monte Carlo method. 

Since the observed qualitative differences in the short-time self- diffusion coefficients of 
charged and uncharged particles are most pronounced for deionized charged suspensions, 
we concentrate here on the case = 0. The system parameters used in our calculations 



are typical for suspensions that have been under experimental study [0, |13|. If not stated 
differently, two-body contributions to HI including terms up to (9(r^^'^) are considered 
together with the leading three-body contribution. Figs. |1| and ^ show our results for 




Figure 1: versus cp for a deionized charge-stabilized suspension with a = 90nm, Z = 200, 

T = 294K, and e = 87.0. Solid line: best fit of the numerical results shows fractional (/>-dependence, 
i.e. Hi = 1 — 2.5941^'^^, with exponent ~ 4/3. Also shown is the dependence of on various 
two-body contributions to the HI. Included terms of the two-body expansion of D^^*^ as indicated 
in the figure. Dashed-dotted line: result for hard spheres according to Eq. @. 

Hi and ifj as functions of (crosses). The corresponding results for hard spheres are also 
included in these figures. Evidently, the effect of HI on Hi and ifj is less pronounced for 
charged suspensions. Moreover, we find a qualitatively different 0-dependence of Hi and 
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Figure 2: Results for obtained for a system with parameters as in Fig. and compared with 
the corresponding result for hard spheres given in Eq. (|^. Best fit of the calculated points (solid 
line) has nearly quadratric (/)-dependence, i.e. = 1 — 1.28(/)^'^^, which extends to surprisingly 
large (j). Further shown is the dependence of on various terms of the two-body series expansion 
of D^^. 

for charged and uncharged particles. Whereas for hard spheres the ^-dependence of if* 
and is linear at small we obtain for charged particles, from a least-square fit of our 
numerical results to the form 1 the following interesting results 

Hi = l-at(f)^-^^, at = 2.59 (7) 
HI = 1 - a,0^•^^ a, = 1.28 (8) 

with exponents close to 4/3 and 2, respectively. Eq. (|^) is valid for (p < 0.05, whereas from 
Fig. 1^ it is seen that Eq. (||) is valid even up to < 0.15. The prefactors at and ar are 
found to be nearly independent of Z for Z > 200. This fact is illustrated in Fig. ^ which 
shows results for Hl{(j)) for various values of Z. Notice that due to the Z-independence 
of Hi and iij, the same Eqs. (0) and (||) are recovered when the accurate, but elaborate, 
Rogers- Young integral equation scheme is used for g{r) instead of the RMSA. 

We will now show that the occurance of exponents close to 4/3 and 2, and the Z- 
independence of at and ar can be understood in terms of a model of effective hard spheres 
(EHS) with density-dependent effective diameter aEus > cr, which accounts for the ex- 
tension of the correlation hole. We can identify (Jehs = ^m, where is the position of 
the principal peak of g{r). It is now crucial to note for deionized suspensions that as 
obtained from the RMSA coincides within 3% with the the average geometrical distance 
f = a{7i / {6(I))Y^^ of two spheres. Thus, we have the scaling relation oc f cx 
Here it is important that Z be chosen large enough that the physical hard core of a par- 
ticle is completely masked by the electrostatic repulsion 0. We now approximate g{r) 
by the rdf gEHsij] (pEHs) of the EHS model, evaluated at the effective volume fraction 
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Figure 3: Hi versus </> for various values of the effective charge number Z as indicated in the 
figure. All system parameters except Z as in Fig. |l]. Notice that becomes nearly independent 
of Z fov Z> 200. 



4>EHS = (t'i.^EHS / (^Y ■ When this approximation for g{r-) is used, and if only the leading 
terms in the series expansions of the two-body moblility functions are retained, we obtain 
the results Hi = 1 — atcj)^^^ and HI = 1 — a^^^ with exponents very close to our numerical 
results. Here 



at = 



ar 



~^'rEHS 
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(9) 



and Gehs{z) is the Laplace transform of xgEHsi^) with x = r/aEHS- Notice that (pEHS, 
and hence at and ar, are independent of and Z (> 200) when cehs is identified as r^. 
To obtain a rough estimate of at and a^, we can further approximate gEHsi^) by 6(x — 1), 
and aEHS by f, giving at = 2.33 and = 0.60. By employing the analytic expression 
for Gehs{z; (pEHs) provided by the Percus-Yevick approximation Jl^, we obtain the values 
at = 3.02 and = 1.12, where the value for in particular is rather close to the numerical 
coefficient in Eq. (H). 

Thus, the EHS model suggests that the scaling relations in Eqs. (0,^ found from our 
numerical calculations are caused mainly by the leading terms in the series expansions of 
the two-body mobility functions. To verify this assertion, we have included in Figs. |] and 
1^ results for if* and ifj obtained by neglecting three-body contributions and by truncating 
the two-body series expansions after various terms of increasing powers in {a/r). These 
figures illustrate our finding that, up to = 0.05, the lowest order contributions to the 
translational and rotational two-body mobilities proportional to and r~^, respectively, 
give by far the most important contributions to Hi and H^. Higher order two-body terms 



and the leading-order three-body term become significant only for > 0.05. For if*, these 
higher-order terms are of the same signature and sum up to increasing deviations in Hi 
from Eq. (0) when is enlarged beyond 0.05. With regard to ifj, however, we observe a 
fortuitous partial cancellation between the three-body contribution and the two-body terms 
of order (9(r~®), which are of opposite sign. As a result, Eq. (^) remains valid even up to 
(j) ^ 0.15. We mention that this cancellation can also be understood in terms of the EHS 
model by reasoning similar to that given above for the leading two-body contribution to 
HI [0]. 

It is further interesting to investigate how and HI are influenced by added electrolyte. 
Our corresponding calculations show a gradual transition from the nonlinear scaling rela- 
tions (0JI) to the expressions (0) and (P) when the amount of added salt is increased 
and when the microstructure changes to a hard-sphere-like structure due to the screening 
of the electrostatic repulsion. 

To summarize, we have calculated the translational and rotational short-time self- 
diffusion coefficient of charged suspensions by incorporating two-body and three-body con- 
tributions to the HI. As a major result, we have found for the first time substantially 
different volume fraction dependencies for charged and uncharged particles. We were also 
able to explain the observed differences in terms of an effective hard-sphere model. We 



mention that recent depolarized DLS experiments [jT3[ on deionized suspensions of opti- 
cally anisotropic particles are in good agreement with our predicted result for H^ in Eq. 
(H). With regard to if*, we are not aware of experimental results that are sufficiently precise 
at low to distinguish the 0^/^-behavior from the essentially linear 0-dependence of hard 
spheres. Finally, we point out that interesting qualitative differences between suspensions 



of charged particles and hard spheres exist also with respect to sedimentation |jT5[ and 
long-time self-diffusion [TH . 
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